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1. Introduction 



In the massless quark limit, Quantum Chromodynamics (QCD) has only one free parameter: 
the strong coupling a s . All strong interaction phenomena should be described in terms of this sin- 
gle input. The measurements of a s at different processes and at different mass scales provide then 
a crucial test of QCD. Obviously, the test should be restricted to those processes where perturbative 
techniques are reliable. Moreover, the same definition of a s should be taken everywhere; the MS 
scheme is usually adopted as the standard convention. Since the running coupling is a function 
of energy, one can either compare the different determinations at the different scales where they 
are measured, checking in this way the predicted momentum dependence of the coupling, or use 
this prediction to bring all measurements to a common reference scale where they are compared. 
Nowadays, the Z mass scale is conventionally chosen for such a comparison. In order to assess 
the significance of the test, it is very important to have a good understanding of the uncertainties 
associated with the different measurements. This is not easy because small non-perturbative effects 
can be present in many observables. In addition, some quantities have been computed to a very 
good perturbative accuracy, next-to-next-to-leading order (NNLO) or even next-to-next-to-next-to- 
leading order (N 3 LO), while others are only known at the leading (LO) or next-to-leading order 
(NLO); the resulting values of a s refer then to different perturbative approximations. The theo- 
retical predictions are also affected by the expected asymptotic (i.e., non-convergent) behaviour 
of the perturbative series in powers of a s . Although this is a common problem of Quantum Field 
Theories, it is probably more severe in QCD because the coupling is rather large (at usual energies). 

The most precise determinations of the strong coupling were reviewed in detail in Refs. [|l|, ||] 
and have been recently updated ^ in the 2012 review of particle physics [Q]. I will heavily use 
the information provided in these references, adding my personal biases and updating the summary 
with the most recent developments. 



2. QCD running and threshold matching 

The scale dependence of the QCD coupling is governed by the renormalization group equation 

da s (pL 2 ) 



P(a s ) = £ A, 



K 



d V n=l 

where the MS /3 function is already known to the fourth order [f| g (£3 = 1.202056903 . . .): 
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The explicit dependence with the number of quark flavours Nf shows the need to properly match 
the different QCD^ effective theories when crossing the corresponding quark thresholds; i.e., the 
value of a s depends on the number of 'active' quarks considered. If one quark flavour is removed 
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Figure 1: Scale dependence of a s at one (dashed), two (continuous, blue) and four (continuous, red) loops 
(left). The dotted line is the 4-loop FOPT approximation in powers of a s (M%) =0.1186. The right plot shows 
the 4-loop error evolution taking a s (M|) = 0.1 186 ± 0.0007; the uncertainty increases at lower energies. 



from the Lagrangian (its mass taken to be heavy enough to decouple), the couplings of the original 
theory with Nf + 1 flavours and the resulting effective one with Nf flavours are related through the 
matching condition 
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where M(jj, 2 ) in the running mass of the heavy quark. The coefficients q„ are known up to four 
loops (Jt<4) @|]: 
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Since c\q = 0, the discontinuity in a s is small for /I = M, but a larger effect is found with other 
choices of the matching point. Physics should of course not depend on this choice and one should 
carefully evaluate the theoretical uncertainties, allowing /» to vary within a reasonable range. 

The 4-loop running coupling allows us to perform a resummation of N 3 LO logarithms, i.e., 
corrections of the form a v (/» 2 )" +4 log" (<2 2 //» 2 ), which are very sizeable when the difference of 
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scales is large. This is illustrated in Fig. [TJ, which shows the scale dependence of a s over a wide 
range of energies, at different levels of approximation. The one-loop j8 -function resummation does 
already a very good job, and the accuracy rapidly improves at higher orders. However, a naive 
fixed-order perturbative expansion (FOPT) in powers of a s (M^) does not approximate well the 
result at low energies, even at four loops, because log [E jMj) is large. 

3. Inclusive observables: cr(e + e — > hadrons) and T(Z — > hadrons) 

The inclusive production of hadrons mediated by y, Z and W ± interactions proceeds through 
the vector V-j = WjY^Wi and axial-vector A^- = YjT^YsYi colour-singlet quark currents = 
u,d,s ...). The theoretical analysis involves the two-point correlation functions 

nj>) = ^^e^(o|r(/JW/x(o)t)|o) = (-^V+^^ngi^+^^n^^ 2 ), 

(3.1) 

where the superscript (0, 1) denotes the angular momentum in the hadronic rest frame. For massless 
quarks, sYlfj^s) = constant (there is a non-perturbative Goldstone-pole contribution to at 

s = 0, which cancels in IT^ 1 '). These correlators are already known to ^(c^ 4 ) [^, [lOt |. In the so- 
called non-singlet contributions (i 7^ j), one quark loop connects the two external currents. If quark 
masses are neglected, these contributions are identical for the vector and axial-vector correlators: 
H(s) = Tlf^jy(s) = YljVj^(s). They are conveniently parametrized through the Euclidean Adler 
function (Q 2 = — q 2 ) 

D(e , )s _ e ^n (e ^^{ 1+£ ,„(^>)"}, 

with Nc = 3 the number of quark colours. The known (n < 4) coefficients take the values [H], [l2|] : 

K y = 1, K 2 = 1.98571 — 0.115295 Nf, K 3 = 18.2427 -4.21585 7V/ + 0.0862069 /v^ , 
K A = 135.792-34.4402/^ + 1.87525^-0.0100928/^. (3.3) 

For neutral currents (i = j) there are additional singlet contributions where each current couples to 
a different quark loop. Owing to the gluon quantum numbers (colour octet and J PC = 1 ), these 
contributions start at &{(X^) and &{a 2 ), respectively, for the vector and axial-vector currents: 

~ £ L « (^)" . ^ - £ £ < (Of)' . (3.4, 

The vector-current coefficients are d\ = -0.41318 and d\ = -5.94225 + 0.191628 Nf [JToJ]. 
The electromagnetic hadronic production in e + e~ annihilation is given by 
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The sum over quark electric charges of different signs strongly suppresses the singlet contribution, 
which has been included as a small correction to the coefficients F„>3. For Nf = 5 flavours, one 
gets Fi = l,F 2 = 1.4092, F 3 = -12.805 and F 4 = -80.434 [0. 

The calculated perturbative series is actually an expansion in powers of CC s (^l 2 ) with coeffi- 
cients containing a polynomial dependence on log(s/fX 2 ). In Eq. (^5|), these logarithms have been 
resummed into the running coupling by taking jj, 2 = s. The physical ratio R e + e - (s) is of course in- 
dependent of the renormalization scale jU, but the truncated series contains a residual /I -dependence 
of &{a^ +1 ), where A = 4 is the last included term, which must be taken into account in the theoret- 
ical uncertainty. Since non-perturbative corrections are suppressed by A 4 /s 2 (the gauge-invariant 
operators contributing to the current correlators have dimensions D > 4), at high energies one can 
perform a N 3 LO determination of a s (s). Unfortunately, the experimental uncertainties are large. 

A much more accurate determination can be obtained from the precise experimental mea- 
surement of the hadronic Z decay rate. The electroweak neutral current j£ = £/(v/V« +afA^) 
involves the vector and axial-vector currents, weighted with the corresponding Z couplings. Since 
cif ■ = 21 f, the singlet axial contributions of the two members of a weak isospin doublet cancel each 
other for equal quark masses; however, the large top mass generates very important singlet axial 
corrections which start at &(a 2 ). The resulting QCD series takes the form 



Rz = 
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(3.6) 



with F x = 1, F 2 = 0.76264, F 3 = -15.490 and F 4 = -68.241 Q10p. Taking properly into account 
the electroweak corrections, the ratio Rz is included in the global fit to electroweak precision data, 
which results in a quite accurate determination of a s (M|). Note that this assumes the validity of the 
electroweak Standard Model with the minimal Higgs mechanism. The PDG quotes the result p|| 



af f - 5) (M 2 ) = a s {M 2 ) 



0.1 197 ±0.0028. 



(3-7) 



4. The hadronic x decay width 



The inclusive hadronic decay of the x lepton provides a very clean way to determine CC S at 



low energies [|13|]. The QCD correlation function of two left-handed W currents receives only 
non-singlet contributions. Restricting the analysis to the dominant Cabibbo-allowed decay width, 

T[x- -> V T + hadrons (S = 0)] 
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Although at low s the integrand cannot be predicted from first principles, the integral itself can 
be calculated systematically by exploiting the analytic properties of the current correlators. They 
are analytic functions of s except along the positive real s-axis, where their imaginaiy parts have 
discontinuities. 7? t ,v + a can then be written as a contour integral in the complex s-plane running 
counter-clockwise around the circle \s\ =m 2 [13-[I"5j]: 
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In the chiral limit (m u j = 0), only the correlator n^^j^s) contributes. Using the Operator Prod- 
uct Expansion (OPE), n( 0+1 )(.s) = Y<dCd/{— s) d I 2 , to evaluate the contour integral, Rt,v+a can be 
expressed as an expansion in powers of 1 /m 2 Jl5|]. The theoretical prediction can be written as 



R%,v+a — N c \V u d\ 2 S EW { 1 + §p + <5np} 



(4.3) 



where Sew = 1.0201 ±0.0003 contains the electroweak radiative corrections [16-|l8|]. The non- 
perturbative contributions are suppressed by six powers of the % mass [15 1 1 and can be extracted 



from the invariant-mass distribution of the final hadrons [19-22]. The presently most complete and 
precise experimental analysis, performed with the ALEPH data, obtains 5np = —0.0059 ±0.0014 



[23]. Quark mass effects [15, 24j 25] amount to a negligible correction smaller than 10 4 . 
The dominant correction (~ 20%) is the perturbative QCD contribution (Nf = 3) 

5p = £^AW(a v ) = a T + 5.20a 2 + 26a T +m a T + --- , 

n=l 

where the contour integrals [|(J 

27ii J\ s \= m \ s \ 71 ) \ m\ m% m% 



(4.4) 



(4.5) 



only depend on a x = a s (m 2 )/7i. The main uncertainty originates in the treatment of higher-order 
corrections []9[ ^3|, ^7]-[36|]. The FOPT expansion in powers of a T , shown in the r.h.s of Eq. (4.4), 
has coefficients much larger than K n because the long running of a s along the circle generates large 
logarithms, log(— s/m 2 ) = ity (0 € [— 71, %]), giving rise to a non-convergent series. These large 
corrections can be resummed to all orders, using the renormalization-group /3-function equation to 
compute exactly (up to unknown j8„>4 contributions) the functions A^ n \a s ). One obtains then a 
well-behaved expansion, known as contour-improved perturbation theory (CIPT) [|^, p7|] , which 
results in smaller values for 5p than the FOPT approach. However, assuming that the known n < 4 
terms of the Adler series are already governed by low-lying infrared renormalons (i.e., already 
sensitive to the asymptotic series regime), it has been argued that CIPT could miss cancelations 



induced by the renormalonic behaviour []28|]; in that case, FOPT could approach faster the Borel 
summation of the full renormalon series. Making optimal conformal mappings in the Borel plane 
and properly implementing the CIPT procedure within the Borel transform, one also obtains nu- 



merical results close to the FOPT value [29, 30]. 

The present experimental value R zy+A = 3.4671 ±0.0084 [||] implies S P = 0.1995 ±0.0033. 
Using CIPT one gets a s (m 2 ) = 0.339 ±0.013, while FOPT would give a s (m 2 ) = 0.318 ±0.014 



[27]. Combining the two results, but keeping conservatively the smallest error, we get 



0.329 ±0.013, 



a, (Ml) 



0.1198±0.0015. 



(4.6) 



Although ag Nf 3 \m 2 ) is significantly larger (16a) than the result obtained from the Z hadronic 



width in Eq. (p.7|), after evolution up to the scale Mz there is an excellent agreement between 



t>ince (l-x) 2 (l+2x) = l-3x 2 +2x 5 , the only non-perturbative contributions to the circle integration originate 
from operators with dimensions D — 6 and 8 (up to tiny logarithmic running corrections). 
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the two measurements. These two determinations provide a beautiful test of the predicted QCD 
running; i.e., a very significant experimental verification of asymptotic freedom. 

At the presently achieved accuracy, a better experimental assessment of 5np would be wel- 
come, which requires a more precise determination of the inclusive hadronic distribution. This 
would also allow for an investigation of uncertainties related with the use of the OPE near the time- 
like axis (duality violations) [23, fw]]; in R r they are heavily suppressed by the presence in (1.2) 



of a double zero ats = m 2 . [15], but they could be more relevant for other moments of the hadronic 
invariant-mass distribution. A recent fit to rescaled OPAL data, with moments chosen to maximize 
duality violations, finds <5np = -0.003 ±0.012 and a s (m\) =0.333 ±0.018 (CIPT+FOPT) [gj, 
in agreement but less precise than the result obtained from the ALEPH invariant-mass distribution, 
Snp = -0.0059 ±0.0014 and a s (m 2 r ) = 0.344 ± 0.009 (CIPT) p3|]. 



5. Lattice determination 

Lattice simulations determine the strong coupling by measuring various short-distance quan- 
tities (non-perturbatively), through a numerical evaluation of the QCD functional integral, and 
comparing the results with the corresponding perturbative expansions in powers of a s , using lattice 
QCD perturbation theory which includes lattice-spacing artifacts. Modern simulations are done 
with 2 + 1 flavours of sea quarks (one tuned to the strange quark and the other two taken with 
masses as small as possible for up and down) and have a NNLO perturbative accuracy. At least one 
dimensionful physical quantity is needed to convert from lattice units to GeV, i.e., to fix the scale 
at which a s is measured. 

The HPQCD collaboration has extracted the coupling, following two different approaches with 
different systematics. They use staggered fermions and the lattice spacing is determined from a 
wide variety of physical quantities. From 22 different simulations of small Wilson loops they obtain 
the value a 4 (M|) = 0.1184 + 0.0006 p2|[. In the second approach, they measure four moments of 
the correlation function of two heavy-quark currents, at 8 different values of the heavy-quark mass 
between m c and m& and 5 different lattice spacings, and obtain a s (M|) = 0.1183± 0.0007 p2|[. An 



independent perturbative analysis, using the results of a previous HPQCD-UKQCD simulation [43] 
already superseded by the new data, finds a slightly larger value Qfc(Af£) = 0.1 192 ±0.0011 Q. 



The PACS-CS collaboration adopts the so-called Schrodinger functional scheme to carry out a 
non-perturbative running of the coupling, from the low energy region used to introduce the physi- 
cal scale to high energies where the matching to the MS scheme is performed. The simulations are 
done with the Iwasaki gauge action and the non-perturbatively improved Wilson-fermion action 
with the clover term. They quote the result a s (Mf ) = 0. 1205 ± 0.0008 ± 0.0005 + [f*|] . The 
Schrodinger functional scheme has been also used by the ALPHA collaboration to study the run- 
ning of a s with Nf = 0, 2 dynamical flavours, but the results with Nf = 4 are not yet complete [|?6|]. 

A numerical simulation of the Adler function, performed with dynamical overlap fermions 
by the JLQCD collaboration, finds a,(M|) = 0.1181 ± 0.0003+°;$^ [0. The uncertainty is 
dominated by the lattice spacing due to the relatively small and coarse lattice used. 

More recently, the ETM collaboration has simulated the ghost-gluon vertex, including a dy- 
namical charm quark (Nf = 2 + 1 + 1), and has studied the running of a s over a wide momentum 
window between 1.7 and 6.8 GeV. Fitting the data with a MOM-like Taylor coupling plus a 1/p 6 
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Figure 2: Lattice determinations of a v (M|). The yellow band shows the average in (5.1). 



non-perturbative term, they obtain a^ Nf 4 \m 2 ) = 0.339 ±0.013, in good agreement with the T 
determination in Eq. (O). They finally quote a s {M%) = 0.1200 ± 0.0014 An analogous 



MOM-like study of the triple-gluon vertex with domain-wall fermions and the Iwasaki gauge ac- 



tion has been presented by RBC/UKQCD at the last Lattice conference [|49j], with the preliminary 
value Ofe(Af§) = 0.1202 ±0.0011 ±0.0002 ±0.0039 ±0.0060, where the last and dominant uncer- 
tainty accounts for the 5% estimated finite- volume error. 

The energy between a static quark and a static antiquark separated by a distance r has been 
also calculated in the lattice, with 2+1 flavours, combining a tree-level improved gauge action 



with a highly-improved staggered quark action [50]. The simulations have been done over a wide 
range of gauge couplings, corresponding to a lattice spacing 1.909/ro < aT x < 6.991/ro, with 
ro = 0.468 ± 0.004 fm. The short-distance part of the static energy can be computed perturbatively 
and it is nowadays known at next-to-next-to-next- to-leading logarithmic (N 3 LL) accuracy, i.e., 
including terms up to order a£ +n log n OC s with n > f^T] — 54 ] (the loga s terms are generated by 



virtual emissions of ultrasoft gluons). Assuming that QCD perturbation theory (after canceling a 
renormalon contribution) is enough to describe the lattice data at distances r < 0.5 ro, a fit to the 
lattice results leads to r A m = 0.70 ± 0.07, or equivalently, ai W/=3) (1.5 GeV) = 0.326 ± 0.0019 



and 0£y(A/§) = 0.1156^q qq22 I55J- Since the present lattice data are not precise enough to profit 
from the N 3 LL calculation, the determination has been obtained at three loop plus leading ultrasoft 
logarithmic resummation accuracy. 

The results of Refs. [4-2, 15, 17, ^8, 55] are largely independent. Following the prescription 
advocated in Ref. []3|], we take a weighted average with a correlated error so that the overall % 2 
equals unity per degree of freedom. This gives 



a s (Ml) = 0.1 186 ±0.0007. 



(5.1) 
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6. Quarkonium 

The bound states of a heavy quark and a heavy antiquark can be rigorously described with 
non-relativistic QCD (NRQCD) techniques, through a combined expansion in powers of a s and 



the heavy-quark velocity v. The strong coupling can then be extracted from the ratio [q6J, |57|] 

r[T(15) -> 7+hadronsl 

R r = Wvf *\1Za T 1 = °- 0245 ± °- 0m ± °- 0013 • ( 6 - 1 

r[T(15) — > hadronsj 

Including the color-octet contributions, which are important in the upper end-point region of the 



photon spectrum, one gets a good description of the photon energy distribution [57]. In order to 



extract a s , a sensible determination of the NRQCD hadronic matrix elements is needed. Using a 



combination of lattice and continuum inputs, the most recent analysis obtains [58] 

«v(M| (15) ) = 0.1841°;™, a s (M\) = 0.119™. (6.2) 

Since the accuracy of this determination is only of NLO in a s (mfy and v 2 , it will not be taken into 
account in our final average. However, the result nicely agrees with the more precise measurements 
performed at other scales. 

7. Particle distribution functions 

The measured scaling violations in particle distribution functions (PDFs) provide precise deter- 
minations of a s , taking advantage of the available deep-inelastic-scattering (DIS) data over a wide 
range of energies and, in particular, the accurate data sets obtained by the HERA experiments. 



A combined analysis of non-singlet structure functions from DIS experiments [p9[], based 
on NNLO and N 3 LO QCD predictions, gave a,(M|) = 0.1 142 ±0.0023, including a theoretical 
error of ±0.0008 This determination neglects possible singlet contributions for x > 0.35, 
where the valence approximation is used. More recent analyses take into account both singlet and 
non-singlet structure functions, together with Drell-Yan and di-muon data, needed for a correct 
description of the sea-quark densities, finding at NNLO the results a s (M%) = 0. 1 134 ± 0.001 1 [|6(J 
and a,(M|) = 0.1158 ±0.0035 Ifj. 

Global PDF analyses include a much broader set of data from fixed-target experiments, HERA 
and the Tevatron. In addition to scaling violations, the determination of the strong coupling ob- 
tained in this way exploits also the dependence on a s of the hard-scattering matrix elements of 
the different processes analyzed. The inclusion of proton collider data allows for a better con- 
trol of the gluon PDF. However, while full NNLO is employed for the structure functions, only 
NLO predictions are available for jet production at hadron colliders. The MSTW group finds 
a v (M|) = 0.1171 ±0.0014 ± 0.002 [|62|], where the last error indicates a conservatively estimated 
theoretical uncertainty. If the Tevatron data is removed from the fit, the central value decreases 
to (X s (Mz) = 0.1104; this DIS-only determination implies a gluon PDF incompatible with the 
Tevatron data, questioning the accuracy of previous DIS fits. In fact, fixing the high-* gluon pa- 
rameters to agree with the global PDF determination, the DIS-only fit gives a much larger value 
« S (M|) = 0.1172 [|2|]. Removing the old (but precise) BCDMS DIS data Q, the result from 
the global PDF fit increases to a s (M|) =0.1181 while the DIS-only fit gives an even larger value 
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a v (M|) = 0.1193 [52]. Pinning down the gluon PDF through Tevatron data seems to be crucial for 
a reliable a s determination, at the expense of some deterioration in the fit quality of BCDMS data. 

More recently, the NNPDF collaboration has performed a global PDF analysis, using neural 
networks as unbiased interpolating functions coupled to a Monte Carlo approach, with the result 
a,(Mf ) = 0.1173 ±0.0011 Q. Including only DIS data in the fit, they find also a smaller central 
value a 4 (M|) = 0.1166, in line with the MSTW results. They observe that at low a s the DIS 
BCDMS data and jet data pull the gluon PDF in opposite directions, so the fit quality can be 
improved using DIS data only in a way which is forbidden when jet data are present. 

As an educated weighted average of all these determinations, the PDG quotes: 

a s {Ml) = 0.1 151 ±0.0022. (7.1) 

The strong coupling has also been determined with inclusive jet measurements at HERA, over 
a wide range of values of the energy scale. The combined result, a s (M§) = 0.1198 ±0.0032 [f|], 
includes a large theoretical uncertainty of ±0.0026. This determination has only a NLO accuracy 
and is not included in the average. 

8. Jet rates and event shapes in <?+e annihilation 

Jet rates have a high sensitivity to a s , which increases with the jet multiplicity n (R n ~ a"~ 2 ). 
Moreover, there are many additional jet observables available, such as a variety of event shapes 
and energy correlations. However, the physical description of jets always involves several scales 
( Emm , px, nib ■■■), making necessary a careful resummation of enhanced logarithmic corrections 
and a good control of non-perturbative power corrections and hadronization effects. Modern studies 
incorporate NNLO corrections [ ]66| , [57| ] and the most recent ones include a matched next-to-leading 



logarithmic (NLL) resummation [68, 69]. Using soft-collinear effective theory (SCET) techniques, 



this resummation has been achieved up to NNLL for jet broadening [70] and N 3 LL for thrust []7l]] 



and heavy-jet mass []72[]. 

A NNLO analysis of the 3-jet rate, using ALEPH data collected at LEP between 91 and 209 
GeV [[73j|, resulted in a, (Mf) = 0.1175 ±0.0025 @. The 5-jet rates at LEP have been also 



analyzed at NLO with the result a s (M|) = 0.1156+q oq3{ [Q. More recently, a re-analysis of the 
3-jet JADE data, at centre-of-mass energies between 14 and 44 GeV, using NNLO predictions plus 
NLL resummation, obtains a s (M|) = 0.1199 ±0.0060 [75], where the main uncertainty originates 



from applying different Monte Carlo models to estimate the transition from partons to hadrons. An 
older study of the 4-jet JADE data, at NLO plus NLL, gave a s (M|) = 0.1 159 ± 0.0028 fiR&. 



Two analyses of six event-shape variables with LEP data, performed at NNLO plus NLL, 



have given the values a s {M£) = 0.1224 ± 0.0039 @ and CC S (M%) = 0.1189 ± 0.0043 [|7g]. In 



both cases, the dominant uncertainty is theoretical. A similar event-shape analysis of JADE data 
gives a s (Mf) =0.1 172 ±0.0051, where hadronization (±0.0035) and perturbative QCD (±0.0030) 
are the largest errors [30|. Including a N 3 LL resummation into the NNLO predictions, a fit to 
the thrust data from ALEPH and OPAL leads to a 4 (M|) = 0.1172 ± 0.0022 0. Using only 
ALEPH data and working also at NNLO plus N 3 LL, Ref. [J72{] finds a similar result from thrust 



and a larger value from the heavy-jet mass distribution, with a final combined result a s (M|) 
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f2 



Dissertori (3j) 
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OPAL (shapes) 
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Abbate (T) 
GLM (T) 
GJL (shapes) 



a s (M z ) 



0.115 



0.120 



0.125 



Figure 3: a s (Mz) from PDFs (left) and e + e jets (right). The yellow bands indicate the PDG []3|] averages. 



0.1193 ±0.0027. Ref. J71| ] includes a ±0.0012 hadronization uncertainty, while hadronization is 
neglected in Ref. [[7^] . 

In the previous analyses hadronization effects are estimated with Monte Carlo generators, 
which were tuned to LEP data. Higher-order perturbative contributions, non-perturbative power 
corrections and hadron-mass effects get then unavoidably mixed, which can have a significant 
impact on the a s extraction. Using SCET techniques, Refs. [ |8~l| ] incorporate explicitly the leading 
non-perturbative power collections, which are also fitted to the data. They analyse the world data 
on thrust distributions at NNLO plus N 3 LL accuracy, including a sophisticated infrared renormalon 
subtraction, obtaining quite low values of the strong coupling. The tail of the thrust distribution 
gives a 5 (M|) = 0.1135 ±0.0011, while cumulant moments using the full thrust distribution lead to 



a s (M|) =0.1 140 ± 0.001 5 [81]. In both cases, the inclusion in the fit of an inverse power correction 



results in a large decrease of the central value (Aa s ~ —0.009), while the total uncertainty has been 
reduced by a factor close to 3 after the renormalon subtraction. The size of subleading power 
corrections, not included in the fit, remains to be investigated. 

A small a s value with a larger uncertainty, a s (M|) = 0.1 131 ^qq^' nas been also obtained 
in another thrust analysis at NNLO plus NNLL [p2||, where hadronization effects are taken into 
account analytically through an effective coupling frozen in the infrared. A previous thrust analysis 
within the same hadronization approach at NNLO plus NLL found a s (M|) = 0.1164^032 [§3|]- 
This approach has also been applied at NNLO to moments of other event-shape variables with the 
result a s (M|) =0.1153 ± 0.0029 @. 



As an educated average of e + e results based on NNLO predictions, the PDG quotes [||] 



a s {M 2 z ) 



0.1 172 ±0.0037. 



(8.1) 



9. Jets at hadron colliders 



The DO collaboration extracted the strong coupling from the pj dependence of the inclusive jet 
cross section in pp collisions at yfs = 1.96 TeV, which is known to NLO, with the result a s (M|) = 



0.1161 lo'oo^g [85]. The uncertainty is dominated by experimental errors (jet energy calibration, pj 
resolution and integrated luminosity). More recently, from a NLO study of jet angular correlations 
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a s (Mz) 



t decay 
Lattice 
PDF 

e + e~ jets 
Z 

Average 



0.114 0.116 0.118 0.120 0.122 
Figure 4: Summary of a s (M|) determinations. 



0.1198 ± 0.0015 
0.1186 ± 0.0007 
0.1151 ± 0.0022 
0.1172 ± 0.0037 
0.1197 ± 0.0028 

0.1186 ± 0.0007 



over a wide range of momentum transfers from 50 to 400 GeV, DO obtains a s (M%) =0.1191 0071' 
finding good agreement with the predicted QCD running up to 400 GeV [j§6|]. 

A first NLO determination based on LHC data was performed with the inclusive jet cross 
section measured by ATLAS at *Js = 7 TeV, testing the running of a s up to jet pj of 600 GeV, 
with the result a,(Af|) = 0.1151 jTo'ooli f@- CMS has made public two preliminary results: the 
measured ratio of the inclusive 3-jet and 2-jet cross sections at y/s = 7 TeV, as a function of the 
average transverse momentum of the two leading jets up to 1.4 TeV, gives a s (M|) = 0.1143 ^5oo67 
[|88[]. while the comparison of the tt production cross section with approximate NNLO calculations 
results in a s {M%) = 0.1178±g;gg^ @. 



10. Summary: the world average value of a s 

Making a world average over the different determinations of a s is a highly non-trivial task 
because systematic uncertainties dominate the most precise measurements. Thus, one relies in 
the more or less conservative attitude adopted to estimate the errors of a given determination. 
Different levels of theoretical accuracy exist for the different observables analyzed and, moreover, 
many theoretical and experimental inputs are highly correlated. Following the procedure suggested 
in Refs. [Jl]-g], we just combine the pre-averages given in the previous sections for each class 
of measurements, which only take into account determinations of at least NNLO accuracy. The 
weighted average of the different input values gives the central value, while the uncertainty is 
adjusted so that the ^ 2 /dof equals unity. One finds in this way 

a s (M|) = 0.1 186 ±0.0007, (10.1) 

which is very close to the value quoted in the 2012 PDG review [Sh. As nicely discussed there, 
the central value remains stable when removing any one of the five inputs included in this average. 
The overall uncertainty is however largely determined by the precise lattice result. 
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